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ABSTRACT 

We use cosmological simulations to identify dark matter subhalo host candidates of the Fornax 
dwarf spheroidal galaxy using the stellar kinematic properties of Fornax. We consider cold 
dark matter (COM), warm dark matter (WDM), and decaying dark matter (DDM) simulations 
for our models of structure formation. The subhalo candidates in CDM typically have smaller 
mass and higher concentrations at z = 0 than the corresponding candidates in WDM and 
DDM. We examine the formation histories of the ~ 100 Fornax candidate subhalos identified 
in CDM simulations and, using approximate luminosity-mass relationships for subhalos, we 
find two of these subhalos that are consistent with both the Fornax luminosity and kinematics. 
These two subhalos have a peak mass over ten times larger than their z = 0 mass. We suggest 
that in CDM the dark matter halo hosting Fornax must have been severely stripped of mass 
and that it had an infall time into the Milky Way of ^ 9 Gyr ago. In WDM, we find that 
candidate subhalos consistent with the properties of Fornax have a similar infall time and a 
similar degree of mass loss, while in DDM we find a later infall time of ~ 3 — 4 Gyr ago 
and significantly less mass loss. We discuss these results in the context of the Fornax star 
formation history, and show that these predicted subhalo infall times can be linked to different 
star formation quenching mechanisms. This emphasizes the links between the properties of 
the dark matter and the mechanisms that drive galaxy evolution. 


1 INTRODUCTION 


The temperature anisotropies in the cosmic microwave back¬ 
ground (CMB) anisotropy spectrum me asured by the Wilki nson 
Microwave Anisotropy Probe (WMAPt dHinshaw et al.ll2013h and 
PLANCK jPlanck Collaboration et alj|20l4 . and observations of 
the large-scale 0.1 Mpc h~^) galaxy clu stering spectrum mea- 
sured by the 2dF Galaxy Redshift Survey |Cole et al.ll2005h and 
Sloan Digital Sky Survey ( ITegmark et al.ll2006ll have shown that 
the large-scale structure for mation is consistent w ith the A cold 
dark matter (CDM) model dFrenk & Whitell2oT3) . The observa¬ 
tion of the CDM component implies physics beyond the stan¬ 
dard model, and many dark matter candidates exist within ex¬ 
tensions t o the standard model of particle physics that behave 
as CDM jjungman et alj Il996h . Though CDM is theoretically 
well-motivated , there are both theoretical (e.g. Abazajian et al.| 
(l 2012 h; Iz ureki (l20l^) and obser vation al dWeinberg et alJ [J2013l: 
iBulbul et al.ll20l4 Boyarsky et al.ll201^ : [Boyarsky et alJl2014 ) in¬ 


terests in considering alternatives. In fact, a broad exploration of 
particle dark matter candidates finds that many viable models be¬ 
have differently than CDM, particularly on small scales, implying 
that observations of dark matter structure on small scales may pro¬ 
vide a unique test of different particle dark matter candidates. 

Aside from the Magellanic Clouds, the eight largest satellite 


galaxies of the Milky Way are dwarf spheroidal (dSph) galaxies 
and the internal kinematics of these dSphs offer one of the best 
prospects for understanding the properties of particle dark mat¬ 
ter from small s c ale as tronomical observations (for a recent re¬ 
view see IWalk^ ( l2013h ). Stellar kinematics unambiguously indi - 
cate that the dSphs are dark matter-dominated dWalker et al.l2007h . 
and their measured potentials have been used to determine whether 
their dark matter profiles are c onsistent with an NFW dens ity pro¬ 
file long predicted by CDM dNavarro et all 1 19961 Il997h . FIow- 
ever, in spite of the high quality data sets, at present it is unclear 
whether the data indicate that the dark matter distributions in dSphs 


are in conflict with the NFW model dWalker & Penamito 


! Amorixco_et_^ 20131) o r are consistent with it ( Jardel et al 
Breddels & Heimj2013l : lstrigari et al.ll2014h . 


2011 


2013 


Improvements in N-body simulations and in hydrodynamic 
simulations of Milky Way-like dark matter halos and their cor¬ 
responding populations of subhalos now provide even more de¬ 
tailed predictions for the dark matter distributions of the dSphs. De¬ 
tailed fitting of the stellar kinematics and photometry to subhalos 
in CDM N-body simulations indicate that the dSphs reside in dark 
matter ha los with maximum c ircular velocity of approximately 20- 
25 km/s dStrigari et al.ll201^ . More general fits to the dSph stel¬ 
lar kinematics in alternative d ark matter model s imulations such 
as warm dark matter (WDM) jLovell et al.ll2014l) . decaying dark 
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matter (DDM)_(|]^n2_e^ay|2014^jand self-interacting dark mat¬ 
ter (SIDM) l lVogelsberger et al~ 1201 2h indic ate that the maximum 
circul ar velocities are larger than in CDM ( iBovlan-Kolchin et alJ 
l201lh . Recent CDM hydrodynamic simulations find that these dark 
matter-only simulations neglect the important effect of baryons, 
which modify the 2 = 0 maximum c i rcular velocities of 

dSphs by about 15% ^ Zolotov_et^J 20121: Brooks_et^J_ ^13l : 


[Brooks & Zoloto'Jl2014l : Sawala et al. 12014 jsawala et al J2OI4I) . 

Identifying subhalos that are consistent with hosting dSphs 
in CDM as well in alternative dark matter scenarios can shed 
light on the cosmological evolution of those subhalos and, per¬ 
haps, the galaxies that they contain. Such identifications may pave 
the way for the development of specific predictions of both CDM 
and alternative models that can serve as true tests of the mod¬ 
els. Exploiting stellar kinematics for this purpose is complemen¬ 
tary to and significantly more robust than using predictions for 
the luminosities of galaxies within subhalos, because these pre¬ 
dicted luminosities are extremely uncertain and rely on extrapolat¬ 
ing phen omenological scaling relations outside of their established 
domain jGarrison-Kimmel et al.ll2014l) . It is additionally comple¬ 
mentary to methods that utili ze measurements of the orbital mo¬ 
tions of Milky Way sat ellites jRocha et al.ll2012l : ISohn et al.ll2013l : 
iKallivavalil et al.ll2013h . 

In this paper we identify dark matter subhalo host candidates 
of the Fornax dSph in CDM, WDM, and DDM N-body simula¬ 
tions by matching to the observed kinematics and photometry. The 
non-CDM based models that we study exhibit DM free-streaming 
effects that are not present in CDM. For WDM models, previ¬ 
ous studies suggest that an equivalent thermal relic mass ~ keV 
generates a tru ncation of the m atter power spectrum on scales 
~ a few Mpc llBode et al.ll200lh . This effect suppresses the for¬ 
mation of small structure below the WD M free-streaming s cale, 
resulting in delayed formation of halos. iLovell et al.l (l2014h use 
galactic zoom-in simulations to show that typical galactic sub¬ 
halos are less concentrated than their CDM counterparts because 
they form at later times. Sterile neutrinos are a canonical WDM 
candidate, and the decays of sterile neutrinos provide a possible 
origin for the detection of an unexplained X-ray line observed 
at 3.55 keV in the Galactic Centre, M31, and galaxy cluster s 
dBulbul et al.l 12014 iBovarskv et all I20l4 iBovarskv et l2014h . 
In DDM models, the DM free-streaming is delayed since the ex¬ 
cess velocity imparted from the DM decay is introduced with a 
lifetime comparable or greater to Hubble time. This can avoid 
the tight l imits placed by h igh-redshift phenomena like Lyman- 


a forest 1 Wang et al. 


_ SS _ 

ture dWang et al.ll2014h . In DDM, at high redshift the structure for¬ 


l2013h and can impact galactic substruc- 


mation is similar to C DM until the age of the Universe is compara¬ 
ble to decay lifetime dWang & Zentnei|[2012h . 

We focus on the stellar kinematics of Fornax because it has 
a high-quality kinematic data sample and the dark matter poten¬ 
tials of subhalos in all of our simulations are well-resolved on the 
scale of the Fornax half-light radius. With the candidate host sub¬ 
halos of Fornax identified, we determine the dynamical properties 
of the host subhalos, such as the present day mass, the peak mass, 
and the maximum circular velocity. With candidate host subhalos 
identified in each simulation, we examine the assembly histories 
of the Fornax host candidates. To assign luminosity to our Fornax 
candidates we use a simple relationship between the stellar mass 
and peak halo mass, and from this we identify Fornax subhalo can¬ 
didates in CDM, WDM, and DDM that are consistent with both 
its kinematics and luminosity. These Fornax candidates allow us to 
predict the infah times and degree of tidal stripping of these subha¬ 
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Figure 1. Top\ DM distribution of subhalos selected from three simula¬ 
tions with different DM models. The navy line represents a subhalo drawn 
from the Aquarius simulation (CDM), the aqua line is from the WDM 
simulation with WDM thermal relic mass = 2.3 keV, and the orange line 
is from the DDM simulations. These three subhalos provide good-fits to 
the Fornax stellar kinematic profile and photometry data. Mass estimations 
for Fornax from previous studies are also shown. The squared data points 
with 1 a error bars sh ow the mass estimation at the 3D half-light radius 
fromlW olf et alj jlOlOl) . and the diamond points with 1 a error bars from 
IWalker et aljl2009l) . Middle'. Slopes of logarithmic mass profiles from these 
three subhalos. The horizontal dash line marks where slope =2, which is 
the theoretical predicted slope for NFW profiles as R —> 0. The non-CDM 
subhalos have slope > 2.0 in the inner region. Bottom'. The best-fit line-of- 
sight velocity dispersion derived using the jeans equation formalism with /? 
= constant model. 


los, and we use these quantities to connect to models for the Fornax 
star formation history in each simulation. 

The outline of the paper is as follows. In § [^ we briefly de¬ 
scribe the properties of simulations used in our analysis. In § [^ 
we review our procedures for identifying Fornax subhalo candi¬ 
dates using its stellar kinematic data and photometry data. In § |4] 
we present our results for the subhalo assembly histories, and dis¬ 
cuss how these quantities can be used to determine the Fornax infall 
time, star formation history, and quenching mechanisms. Lastly we 
draw our conclusions in §[^ 


2 SIMULATIONS 

In this section we briefly describe the cosmological simulations that 
we utilize. For more details we refer to the original simulation pa- 
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Figure 2. The 68% and 95% contour region for the stellar density slope at r=600 pc versus velocity anisotropy /3. From left to right, the contours are drawn 
from the MCMC results of the Jeans equation fitting for the subhalos shown in Figure[T](left: CDM, middle: WDM, right: DDM). 


pers (CDM : Sorineel et alJ 

2008h. WDM 

and DDM :IWang et alj (12014 



For the CDM model, we utilize the Aquarius simulations, 
which are si x realizations (Aqua rius A to F) of galactic zoom-in 
simulations l lSpringel et al.ll20 08h. These simu lations are generated 
using the Gadget code l lSpringel et alj|2008h and use cosmologi¬ 
cal parameters consistent with the one-year and five-year Wilkinson 
Microwave Anisotropy Probe (WMAP) data: Hq = 73 km/s/Mpc, 
Dm = 0.25, Da = 0.75, as = 0.9, and = 1. We adopt the 
level-2 resolution simulations as our main sample (with Plummer 
gravitational softening length e = 65.8 pc and particle mass rrip 
= 1.399X IO^’M©—6.447X lO^M©), and utilize the highest level- 
1 resolution for Aquarius A (with e = 20.5 pc and particle mass 
trip = 1.712X 10^Mq) to perform resolution tests to understand the 
effects of the force softening scale on our analysis. The details of 
these tests are described in the Appendix §|^ The properties of the 
Aquarius simulations are shown in Tablej^ 


For the WDM models, we use the simulatio ns described 
Lovell et H] (l2014ll . To account for WDM physics, ILovell et al] 


(l2014h re-simulate the Aquarius A halo using initial condition 
wave amplitu des that ar e resca led with thermal relic WDM power 
spectra from IViel et al.l l l2005ll . We adopt their “high resolution” 
suite that corresponds to level-2 in the original Aquarius nota¬ 
tion with rrip = 1.55 xIO'^Mq and e = 68.1 pc. This suite in¬ 
cludes two WDM simulations with equivalent thermal relic masses 
of 2.3 keV and 1.6 keV, and their CDM counterpart simulations. 
The 2.3 keV thermal relic is a good approximation to the mat¬ 
ter power spectrum of a 7 keV sterile neutrino that is resonantly 
produced in a lepton asymmetry Le ~20, which translates to a 
transfer function warmer t han that infer red from the 3.55 keV line 
jVenumadhav et alj|2015l : Il7ovellll201^ . We take 1.6 keV model 
as a case that is ruled out by current Lyman-a forest limits (e.g. 
dViel et al.ll2013h j and 2.3 keV as the model likely allowed by such 
limits for comparison. The cosmological parameters are derived 
from WMAP7: Ho = 70.4 km/s/Mpc, = 0.272, Da = 0.728, as 
= 0.81, and ris = 0.967. Note that this is different than the original 
Aquarius simulations for which WMAPl cosmology was imple¬ 
mented. We describe the effects of the different cosmology on our 
results in Appendix §|^ The self-bound halos in Aquarius simula¬ 
tions and WDM simulations w ere identified using the SUB FIND 
algorithm iSpringel et alj200i1) . 


For our DDM model, we use th e set of simulation s that im¬ 
plement late-time DDM physics from lWang et alj ( l2014h . In these 
DDM models, a dark matter particle of mass M decays into a less 
massive daughter particle of mass m = (1 — f)M with / ^ 1 


and a significantly lighter, relativistic particle, with a lifetime F”^, 
where F is the decay rate. The stable daughter particle acquires 
a recoil kick velocity, 14, the magnitude of which depends upon 
the mass splitting between the decaying particle and the daugh¬ 
ter particle. The DDM simul ations are genera ted using a modi¬ 
fied version of GADGET-3 jPeter et al. 1 20101) . Here we consider 
the case with decay lifetime F“^= 10 Gyr and kick velocity 14 = 
20 km/s. This model has been show to ha ve interesting imp lica- 
tions for dark matter small-scale structure (jM^ng^et^y 2014h an d 
is allowed by current Lyman-a forest limits jWanget al.ll2013h . 
The cosmology used is based on WMAP7 results with 770=71 
km/s/Mpc, Dm=0.266, Da=0.734, <78=0.801, and ns=0.963. We 
use Amiga Halo Finder (AHF) jKnollmann & Knebel 12001^ 
for halo findi ng and the merger tree is constructed using CONSIS¬ 
TENT Trees feehroozi et alj2013ll . 

We define the virial mass (M 2006 ) of the galactic halo as the 
mass enclosed within the region of 200 times the b ackground for 
all sim ulations, which corresponds to the M 50 in ISpringel et al.l 
( l2008h. We adopt this definition for the same reason quoted 
in ISpringel et alj (l2008h . namely that it yields the largest radius 
among other conventional halo definition and hence the largest 
number of substructures. 


3 FITTING SUBHALOS TO STELLAR KINEMATICS 


In this section we discuss our method for fitting subhalos in simu¬ 
lations to the stellar kinematic and photometric data of Fornax. For 
our theoretical analysis we use the spherical jeans equation, and 
allow for a constant but non-zero anisotropic stellar velocity dis¬ 
persion. We utilize this theory in order to efficiently identify a large 
sample of Fornax candidates. For co mparison studies with simi¬ 
lar motivations have been und ertaken jBovlan-Kolchin et al.ll201ll : 
iGmison-Kimrnel et al.l l2014h that have use d the mass estimator 
of lWalker et al.l ( l2Q09h and iwolf et al.l ( l201(]h to estimate maximum 
circular velocities for the dSphs. As we show below our results are 
in good agreement with these previous results. Since we start at the 
level of the jeans equation, we are able to identify preferred orbital 
structure of the stars given the underlying subhalo potentials. 

We assume that the potential is spherically symmetric, 
dispersion-supported, and in dynamical equilibrium, so that we can 
derive the stellar line-of-si ght velocity dispersion pr ofile as a func¬ 
tion of projected radius R dBinnev & Mamonlll982h : 
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Simulations 

Particle mass rrip 

Force Softening e 

M2006 

r2006 

Dark Matter Properties 


[Mq] 

[pc] 

[Mq] 

[kpc] 


Aq-Al 

1.712x10^ 

20.5 

2.52x10*2 

433.5 

CDM 

Aq-A2 

1.370xl0'‘ 

65.8 

2.52x10*2 

433.5 

CDM 

Aq-B2 

6.447 xlO^ 

65.8 

1.05x10*2 

323.1 

CDM 

Aq-C2 

1.399x10-* 

65.8 

2.25x10*2 

417.1 

CDM 

Aq-D2 

1.397x10-* 

65.8 

2.52x10*2 

433.2 

CDM 

Aq-E2 

9.593x103 

65.8 

1.55x10*2 

368.3 

CDM 

Aq-F2 

6.776x103 

65.8 

1.52x10*2 

365.9 

CDM 

Aq-A2 w7 

1.545x10* 

68.2 

2.53x10*2 

432.1 

CDM 

Aq-A2-mi.6 

1.545x10* 

68.2 

2.49x10*2 

429.9 

WDM (mwDM=1.6keV) 

Aq-A2-m2,3 

1.545x10* 

68.2 

2.52x10*2 

431.4 

WDM (mwDM=2.3keV) 

Z13-CDM 

2.40x10* 

72.0 

1.31x10*2 

335.2 

CDM 

Z13-tl0-v20 

2.40x10* 

72.0 

1.16x10*2 

336.2 

DDM {r -*=10 Gyr, ^^.= 20.0 km/s) 


Table 1. Parameters of simulations. The mass of the galactic halo M 2006 is defined as the mass enclosed within the region of 200 times of the background 
density. 


with 


h{R) = 2 



\/r'^ — 


( 2 ) 


Here p*(r) is the 3D stellar density profile and I*{R) is its 
2D projection. The velocity anisotropy parameter is /3(r) = 1- 
cTj (r) /a'^{r), where ar{r) is the radial velocity dispersion of stars 
and (Jt{r) is the tangent ial velocity dispersion. The se quantities sat¬ 
isfy the Jeans equation dsinnev & Tremain3l2008ll : 


d{p*(yl) , \ 2 , sGM{<r) 

r -:-h 2p{r)p,a^ = —p*[r) - 


(3) 


dr . • - - 

For the case of constant non-zero /3, the solution of Eq. has a 
simple form: 


p*a^(r) 


1 r p.{r')GM{<r')r'^^ ^ , 

/ - —2 - 

J r 


(4) 


To model the 3D stellar density profile p, (r) as a function of 
radius we use this general form: 


p»(r) oc 


1 

a;“(l -I- a;*’)(=-“)/<>’ 


(5) 


where x=r/ro and a, b, c are free parameters that capture the stel¬ 
lar distribution slopes over different radii. A density profile of this 
form has been f ound to adequately describe the photometry of the 
classical dSphs dStrigari et al.ll201ol) . Though the conversion of 3D 
to 2D profile is not a one-to-one relation, implying that different 
choice of parameters for the 3D profiles can provide similar 2D 
projected profiles, the present photometry data do give good con¬ 
straints on the Fornax stellar distributions with profiles of the form 
Eq.[5] We fix the normalization of stellar mass by assuming mass- 
to-light ration M,/L = 1 a nd adopt the Fornax V-band luminosity 
value of L„=1.7x10^Mq ( lMcConnachiell2012h . 

With the above equations we utilize the following algorithm to 
model the Fornax line-of-sight velocity dispersion given the subha¬ 
los in each of our simulations. We begin by finding subhalos in our 
galactic zoom-in simulations using halo finders, and for each sub¬ 
halo determine the dark matter distribution and thus the potential 
as a function of radius. Given the potential of each subhalo, the 


stellar distribution described by the parameters a, b, c, ro, and ve¬ 
locity anisotropy by /3, we solve the jeans equation to determine the 
line-of-sight velocity dispersion. We then fit to the photometric and 
kinematic data by marginalizing over these parameters a, b, c, ro, 
and f} via a Markov Chain Monte Carlo (MCMC) method to deter¬ 
mine the best fitting value for the model parameters. 

For our fits to the Fornax stellar kinematics we define 

[ai - aios{Ri)f 

Xa - 2^ -3-> (6) 

where Nuns is the number of annuli, Ri is the mean value of the 
projected radius of stars in the annulus, and aioaiRi) is the 
derived velocity dispersion for a given subhalo in the annulus. 
The line-of-sight velocity dispersion from the binned data is Oi and 
the corresponding error is in each annulus. 

The kinematic datasets that we use consist of line -of-sight stel¬ 
lar velocities from the samples of lWalker et alj(l2009h . We consider 
only stars with > 90% probability of membership, which gives us 
a sample of 2409 Fornax member stars. We bin the velocity data 
in circular annuli and derive the mean line-of-sight velocity in each 
annulus as function of Ri. We calculate the line-of-sight velocity 
dispersion d and their error e in e ach annulus following the method 
described in IStrigari et al]( l2010h . 

For our fits to the Fornax photometry we define 

^,2 ^[h{Ri)-h[Ri)]^ 

xiR-y -^2-> G) 

i—1 ’ 

where Nm is the number of radial bins, Ri is the radius of the 
data point, and is the derived 2D stellar surface density 

in the radius bin and ef ^ ^ is the unce rtainty. For our photo¬ 
metric data we use the measurements from IColeman et ^ ( l2005h . 
Note that in our analysis we assume that Fornax contains a single 
population of stars. 

For our MCMC, we adopt uniform priors over the following 
range: 0 ^ a 2, 0.5 $jb^5, 4^c^8, 0.5 ro ^ 2, -1.0 
^ /3 ^ 1.0. The MCMC calculation is performed using a modified 
version of the publicly available code COSMOMC dLewis & Bridla 
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l2002h as an MCMC engine with our own likelihood functions. We 
adopt a simple form of the likelihood function as the summation 
of the from both the stellar kinematic data and photometry over 
their degree of freedom: Xtot/d.o.f. = xl/d.o.f. + Xm/d.o.f., 
where fa and fiR are calculated using Eq. |^and Eq. |7]with de¬ 
gree of freedom (d.o.f.) equal to number of data bins minus the 
number of free parameters plus one. We have 15 radial bins in the 
velocity dispersion data and 19 radial bins in the photometry light 
prohle data. We then select the subhalos that satisfy the condition 
of Xtot/d-o.f. ^ 3.0 as good-fits to the Fornax kinematic and pho¬ 
tometry data. We choose these criteria in order to obtain a conser¬ 
vatively large subhalo sample of candidates that can host Fornax. 


4 RESULTS 

We now present the results of our analysis. We begin by presenting 
the properties of the dark matter subhalos in the different simula¬ 
tions that are candidate Fornax hosts using stellar kinematics alone. 
We then combine with a simple model for subhalo luminosities to 
identify subhalos that match both the luminosity and kinematics of 
Fornax, and present results for the assembly history of these sub¬ 
halos. 


4.1 Properties of Fornax Candidate Host Subhalos at z=0 

In this subsection we discuss the properties of Fornax candidate 
host subhalos at 2 : = 0. We begin by discussing the density prohles. 
for which we compare the central densities of the subhalos from the 
different simulations. We then discuss the total subhalo masses and 
maximum circular velocities of the host subhalo candidates. 


4.1.1 Subhalo density profiles and correlations 

We begin by examining the density prohles of the subhalos that 
provide good hts to the photometry and line-of-sight velocity dis¬ 
persion prohles. Figure [T] shows the line-of-sight velocity disper¬ 
sions for three different representative subhalos, one from each of 
the CDM, WDM, and DDM simulation. For CDM, here we choose 
the subhalo that has the minimum value of Xtot/d.o.f., while for 
DDM and WDM we choose the subhalos that have Xtot/d.o.f. 
< 3 and have the shallowest central density prohles. For WDM 
(DDM), the Xtot/d.o.f. is specihcahy 1.67 (0.85), and the p-value 
for the velocity dispersion ht is 0.4 (0.84). For the CDM subhalo, 
the Xtot/d.o.f. is 0.93 while the p-value for the velocity dispersion 
ht is 0.8. Thus in all three cases the hts are good from a statistical 
standpoint. 

Figure [T] also shows both the mass distribution and the log- 
slopes as a function of radius of the selected subhalos. The CDM 
subhalo is the most centrally-concentrated at radii f 500 pc, gen¬ 
erating a central slope of AlogioM/AlogioR ~ 2, which is con¬ 
sistent with the predicted value from an NFW prohle. The DDM 
subhalo has a central slope of ~ 2.6, while the WDM subhalo 
has a central slope of ~ 2.5. Thus the WDM subhalos have “soft” 
cores, which are expected since WDM subhalos are less concen¬ 
trated than their CDM c ounterparts because they form at a later 
time dLovell et al.ll2012h . The DDM subhalo has a slightly shal¬ 
lower core-like feature than in W DM because the effects of dark 
matter decay are scale-dependent dWang et al.l2014l) . 

As Figure[T]indicates, the predicted line-of-sight velocity dis¬ 
persions are very similar, even though the best htting values of 


Table 2. Average subhalo properties of Fornax kinematic good-fits 


Simulations 

Msub 

^peak 

^max 


{Mq\ 

[Mq] 

[km/s] 

Aq-A2-w7 

3.72±1.91xl0® 

I.55±2.01xl0® 

17.70±1.02 

Aq-A2-m2.3 

1.30±1.32xl0^ 

4.23±2.09xl0® 

20.31 ±2.26 

Aq-A2-mi.6 

3.52±3.15x10® 

5.05±2.67xl0® 

23.22±4.57 

Z13-CDM 

4.50±2.00xl0® 

8.77±2.94xl0* 

18.26±0.71 

Z13-tl0-v20 

2.11±1.78xl0® 

4.87±5.10x10® 

25.03±3.18 

Aq-A2 

3.46±1.27xl0* 

I.60±2.14xl0® 

17.69±0.80 

Aq-B2 

2.34± 1.06x10* 

1.04± 1.76x10® 

17.78±0.98 

Aq-C2 

2.60±1.10x10* 

0.91±1.25xl0® 

17.61±0.44 

Aq-D2 

3.40±2.30xl0* 

8.28±4.23xl0* 

17.81±1.13 

Aq-E2 

2.84±1.23xl0* 

6.62±3.79xl0* 

17.60±0.97 

Aq-F2 

3.58±2.43xl0* 

1.00±0.59xl0® 

17.98±1.04 

Aq Average 

3.07±1.77xl0* 

1.02± 1.29x10® 

17.76±0.94 

Aq Max 

1.18x10® 

8.83x10® 

21.71 

Aq Min 

1.09x10* 

3.08x10* 

15.79 


a, b, c, To and are different in each case. This highlights an ex¬ 
plicit degeneracy between the dark matter potential, the veloc¬ 
ity anisotropy, and the stellar density profile. At the level of a 
jeans-based analysis, the degeneracy between the dark matter po¬ 
tential and the velocity anisotropy has been well-studied and has 
been known to preclude d etermination of the dark matter profile 
shape dStrigari et al.l200% . 

The degeneracy between the stellar density profile and both 
the dark matter potential and the velocity anisotrop y is less well- 
understood (see however lstrigari et alj(l2010l . l2014l) ). so it is inter¬ 
esting to examine this further within the context of a Jeans-based 
analysis. Figure [^highlights this degeneracy for the case of a sin¬ 
gle subhalo by showing the correlation between the stellar density 
slope at r = 600 pc versus anisotropy parameter. For ah models, 
within 68% contour region j3 spans a wide range and is degenerate 
with the stellar density concentration. The specific allowed regions 
change in each of the models, because they predict different DM 
density slopes. 

4.1.2 Subhalo masses and circular velocities 

Though the density profiles of the CDM, WDM, and DDM subha¬ 
los are indistinguishable from the photometric and the kinematic 
data, the subhalos are distinct when considering more global prop¬ 
erties. In Figure 13 we show the circular velocity curves, Vdrc = 
\JGM{< r)/r, of subhalos that provide good fits to the Fornax 
kinematic data and photometry data. Note that for each curve the 
stellar mass contribution has been included, so that at a radius of 
0.9 kpc, the enclosed mass contribution from stars is about 10% of 
the DM mass. The circular velocity curves in the WDM and DDM 
models are shallower in the center than the CDM subhalos, which 
is consistent with the discussion above on the individual best-fitting 
density profiles. Further, the WDM and DDM subhalos have larger 
ma ximum value, Vmax, which is consis tent with previous analy¬ 
ses ( iLovell et alj2012l : IWang et alj2014h . 
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Figure 3. Circular velocity curves of subhalos that are good fits to the Fornax stellar kinematic and photometry data. Notice that for each curve the stellai' 
mass contribution has been included. The dark gray shaded areas include 68% of the subhalo curves from Aquarius simulation A to F, and the light gray areas 
are for 95 % of the subhalos. The solid black lines are the average of all the Aquaiius CDM fits. From left to right panel, the colored lines sho\v the subhalos 
from the WDM simulation with 2.3 keV thermal relic mass (left panel, aqua), the WDM simulation with 1.6 keV thermal relic mass (middle panel, green), 
and the DDM simulation with lifetime =10 Gyr and kick velocity = 20 km/s (right panel, orange) 
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Figure 4. Histograms of subhalo properties for the Fornax good-fits. Here we show two subhalo present-time properties: maximum circular velocity (Vmax, 
left panels), and subhalo mass right panels). In the upper two panels the navy histograms show the Aquarius A results, the aqua ones are for the WDM 

2.3 keV simulation, and the green ones are for the WDM 1.6 keV simulation. In the lower two panels the dark red histograms show the Z13 CDM results, and 
the orange ones are for the DDM simulation. We note that for the same galactic halo realization, the alternative DM scenarios presented here (WDM, DDM) 
generate much less substrucure than their CDM couterpaifs (see text in 1( 14.1.21 . 



Figure|4] shows the distrihution of Vmax and the total mass of 
the suhhalo, Msub, for the CDM suhhalo host candidates that are 
good fits to Fornax. Flere we show only results from the CDM sim¬ 
ulations Aq-A2-w7 and Z13, because from these simulations we 
are able to identify the suhhalo counterparts from the WDM and 
DDM simulations, respectively. From Figure|4]and also Table[2 we 
show that CDM host suhhalo candidates have 10®M© ^ ^ 

10®Mq, 16 km/s ^ Vmax 22 km/s. Table|^shows the average 
Msub and Vmax of Fornax candidates from all our simulations. As 
is indicated, Aq-A2-w7 and Z13 provide a good representation of 
our entire sample of CDM simulations. The Vmax range of candi¬ 


dates from six Aquarius simulations (see Table have an average 
~17 km/s and an associated small variance, which agree well with 
results from previous study dBovlan-Kolchin et al.ll201 j) . 

For comparison Figure |4] shows the WDM and DDM subha¬ 
los that match the Fornax kinematics and photometry. In the WDM 
simulation with 2.3 keV WDM mass we identify seven host sub¬ 
halo candidates with Msub ~ 0.3 —3 x 10® Mq and Vmax ~ 18-23 
km/s. In the WDM simulation with 1.6 keV particle mass, we iden¬ 
tify six host suhhalo candidates with Msub ~ 0.2 —8 x 10 ® Mq and 
Vmax ~ 17-28 km/s. In the DDM simulation we identify three host 
suhhalo candidates with Msub ~ 0.8 — 4 x 10® M© and Vmax ~ 
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Mpeak [M®] Mp„k [Me] Mp„k [Me] 

Figure 5. The Mp^ak distribution of the Fornax good-fit subhalos. The color notation in the left panel is the same as previous figures, and the right panel 
shows the distribution from all six Aquarius level-2 simulations (A-F). We note that there is only one galactic halo realization in the left and middle panel, 
while there are six in the left panel. The y ellow shaded areas ind icate the region where the luminosity of these objects are consistent with the Fornax luminosity 
using abundance matching methods from lBehroozi et ^ j2Q13h . 


22-28 km/s. These results indicate that from our Jeans-based mod¬ 
eling WDM and DDM predict a population of massive Fornax host 
subhalos that have shallower central dark matter density profiles. 
We note that alternative DM scenarios always predict much fewer 
Fornax candidates than their CDM counterpart simulations, sim¬ 
ply because the number of subhalos are reduced significantly due 
to di fferent DM properti es. For example for the WDM case. Fig. 
11 in lLovell et alj ( |2014[) shows that the numbers of subhalos with 
Vmax ^ 15 km/s are 14 (1.6 keV), 28 (2.3 keV) and 120 (CDM). 
For the DDM case the subhalo number with Vmax -5 15 km/s is 16 
and 56 for the CDM counterpart. The subhalo numbers in different 
CDM halos reflect mainly the differences in halo mass and merger 
history. 


4.2 Predictions for Luminosity and Effects of Reionization 

To this point we have not used the luminosity of Fornax as a con¬ 
straint, other than to consider its impact on the stellar kinemat¬ 
ics and photometry. In doing so we have implicitly assumed that 
all of the subhalos are suitable hosts of a galaxy with the present 
day luminosity of Fornax. We have made this assumption in or¬ 
der to analyze a large statistical sample of well-resolved Fornax 
host candidates in N-body simulations that do not account for 
the effects of baryons on the subhalo evolution. High resolution 
Local Group simulations with baryons have been recently under¬ 
taken dSawala et al.l2014| : ISawala et al.l201^ . though at present the 
statistical sample of subhalos is smaller in these simulations than 
the corresponding sample from N-body simulations. For this rea¬ 
son, we use semi-analytic models for the luminosities of our Fornax 
subhalo host candidates, and leave to a future study the analysis of 
the simulations that include the baryons. 

Our subhalo luminosities are motivated by the results of sev - 
eral studies. First, we consider the results of ISawala et al.l ]2014h . 
who utilize hydrodynamic simulations to predict the luminosity of 
subhalos before they fall into the Milky Way halo and z — 0. 
For subhalos with the present day luminosity of Fornax, Sawala 
et al. find that subhalos with the luminosity of Fornax have a to¬ 
tal mass in the range ~ 2 — 10 x 10® Mq, and a total peak 
halo mass in the range ~ 7 — 20 x 10® Mq. The high end of 
this mass range is consistent with previous semi-analytic mod¬ 


els iCooper et alj|201ol) and with the predicti ons from the abun¬ 
dance matching method jBehroozi et al .11201^ . For dSphs fainter 
than Fornax, the predicted range of subhalo masses for a fixed lumi¬ 
nosity is larger, and also there are stronger deviations from the p re- 
dictions of abundance matching and that of I Sawala et alj ( 120141) . A 
well-known caveat, however, is that the abundance matching tech¬ 
nique is not calibrated at both sub-galactic scales and in non-CDM 
cosmologies. Therefore results in these regimes depend on extrap¬ 
olating the existing models and assuming the average galaxy for¬ 
mation history holds for non-CDM based models. 

In order to determine the relevant peak subhalo mass range 
for Fornax, we first take the lower bound on the stellar mass of For¬ 
nax to be Mff lxIO^M q. The lower bound is derived from the 
1 a lower bound on the Fornax V-band luminosity of Lv =1.7 
±g;t X wMq, and assuming a stellar mass-to-light ratio M*/L = 
0.8 (M* ~ 1.3 X 10^ X 0.8Mq for the lower bound). For this 
adopted stellar mass, we then select Fornax candidates by their 
peak mass that predict the sam e stellar mass range us ing the abun¬ 
dance matching description in iBehroozi et alJ ) l2013h . In Figure |5] 
we show the distribution of peak subhalo masses for these Fornax 
candidates. These results indicate that while subhalos that are con¬ 
sistent with the Fornax kinematics have subhalo peak masses in 
the range Ms„(, ~ 10®-10 ^®Mq, only a handful at the very mas¬ 
sive end predict the right luminosity. For example from the right 
panel in Figurel^we find that 2 out of 124 candidates from the six 
Aquarius simulations are consistent with the Fornax luminosity. In¬ 
terestingly, these two candidates come from the Aquarius A and B 
simulations, which are, respectively, the most massive and the least 
massive halos among all six Aquarius simulations. This would in¬ 
dicate that the likelihood of generating a Fornax candidate may be 
only weakly correlated with the galactic halo mass. 

For comparison, in Figure we show the Fornax candidates 
from the WDM and DDM simulations, as well as their CDM coun¬ 
terparts. Though the sample of host subhalos is smaller than in the 
case of CDM, there is a higher probability for WDM and DDM 
Fornax candidates to reside in more massive subhalos. Therefore, 
the WDM and DDM candidates are more likely to match the bright 
luminosity of Fornax than the CDM candidates. In each case, we 
find 1-2 candidates that match the kinematics and luminosity, out of 
3-7 candidates that just match the kinematics. Note again the caveat 
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M,=6 [Mo] M,=e [M®] 

Figure 6. The distribution of the progenitor mass at z=6 for the Fornax candidates. The color scheme in the left panel is the same as Figure|4]showing DDM 
and Z13 simulation candidates, and right panel shows all the candidates from all six Aquarius halos. We note that there is only one galactic halo realization 
in the left panel, while there are six in the left panel. The filled histograms indicate the subhalo s with luminosity tha t match Fornax at z=0. The shaded areas 
indicate the progenitor mass range with baryon content subject to UV background suppression JOkamoto et alj2008l) . 


that the relationship between the stellar mass and subhalo progeni¬ 
tor mass may differ from this relation which is derived from CDM 
simulations. 

With candidate host subhalos identified we are now in posi¬ 
tion to study their evolution and the masses of their progenitors 
near the redshift of reionization. The effect of reionization on sup¬ 
pressing star formation in small halos shows that halos below a few 
10® Mq likely have no stars due to the UV-background suppres¬ 
sion. This implies that any subhalo with a progenit or less massive 
would be unlikely to ho st a visible galaxy today dOkamoto et al.l 
120081: ISawala etal.ll20l4 . 

In Figure we show the subhalo progenitor mass at z=6 for 
the Fo rnax candidates. According to the criteria of lOkamoto et al.l 
ll2008l) . more than half of the Fornax candidate in Aquarius would 
not host any stars. The subhalos that match both the Fornax kine¬ 
matic and photometry data as well as the luminosity are on the 
massive end of the distribution. 

In the left panel of Figure we show the case for the DDM 
simulation and its CDM counterpart. We do not include WDM sim¬ 
ulations here because their structure formation and galaxy forma¬ 
tion at high redshift are expected to be significantly different than 
in CDM. Also for the WDM models that we considered here, ma¬ 
jority of the subhalo progenitors are not yet formed or identified 
by halo finders at z=6. However, for the DDM model with a de¬ 
cay lifetime of 10 Gyr, it is safe to assume a similar reionization 
history as CDM. Since at 2 = 6 less than about 10% of the dark 
matter has decayed with a small recoil kick velocity of 20 km/s, 
the structure growth of both the dark matter and baryon compo¬ 
nents should be similar to CDM. Again the majority of the CDM 
subhalos should be UV-suppressed, and all three DDM subhalos 
are above the threshold, with one object that matches the Fornax 
luminosity. 


4.3 Mass Assembly Histories 

We now move on to discuss the mass assembly histories of our 
Fornax candidates. The left panel of Figure |3 shows the subhalo 
mass as a function of redshift for all 124 of our Fornax candidates 


from all six Aquarius halos. The middle panel shows the corre¬ 
sponding candidates from the WDM model with 2.3 keV mass and 
the right panel shows the candidates from DDM simulations. From 
the discussion above, we identify the subhalos that match the For¬ 
nax luminosity, and also determine an approximate infall time for 
those subhalos that match both the Fornax luminosity and kinemat¬ 
ics. From Aquarius simulations, we can see that the candidates that 
match the kinematics and luminosity, which are shown in blue lines 
rather than the gray lines for candidates only fit stellar kinematics, 
also have the largest ratio of yi{Zpeak) /M {z = 0), with typical 
peak values for this ratio of ~ 30—60 at 2 ~2. For WDM, subha¬ 
los that match the kinematics and luminosity also show a large peak 
value at 2 ~ 1.5. On the other hand the corresponding DDM sub¬ 
halo shows a relatively low peak value of yi{Zpeak)/M[z = 0) ~ 
3 at 2 < 1. 

From the results presented in Figure |7] we can identify three 
important trends and predictions regarding tidal stripping, infall 
times, and star formation histories. 


4.3.1 Tidal Stripping 

Using the above results in CDM we have a clear prediction that For¬ 
nax has lost a significant amount of dark matter mass after it falls 
into the Milky Way. Does this translate into a more specific obser¬ 
vational signature? One possibility is that this tidal stripping of dark 
matter is manifest in the tidal stripping of the stars. However, since 
our models only directly include dark matter, it is not clear whether 
we should expect stellar tidal tails to be observable. It is possible 
that the stellar components are embedded deeply in the center of 
the subhalo potential, so that the dark matter is mostly stripp ed and 
the stars left unstripped. For example, IWatson et al.l ( 1201 21) shows 
that if subhalos experience substantial dark matter mass loss before 
mass is lost within the galaxy, this explains how satellite galaxies 
lose stellar mass and contribute to “intrahalo light” (IHL). 

For comparison the Fornax candidates in WDM, shown in the 
middle panel of Figure |3 are similar to CDM in that both require 
large Mpeak/M{z = 0). In DDM, as shown in the right panel 
of Figure [ 7 ] the ratio is much smaller, Mpeak /M (2 = 0) ~ 3. 
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Figure 7. Mass assemb ly histories for the Fo rnax good-fit subhalos. The gray lines show the fits with luminosity dimmer than Fornax using abundance 
matching technique from iBehroozi et ^ i2013h . and the blue lines show those with luminosity that matches the Fornax luminosity. The vertical dotted lines 
show the infall time of the blue curves. In the left panel the dashed black lines include 68% of the Aquarius simulation fits, and the black solid lines show the 
mean value. 


These results hint at the interesting possibility that hosts of Fornax 
have undergone different degrees of tidal stripping in different dark 
matter models. Thus with future larger statistical samples of non- 
CDM models, as well as the inclusion of baryons in the simulation, 
the degree of tidal stripping of Milky Way-satellites may provide a 
new test of dark matter models. 


4.3.2 Inf all Time 

From the peaks in the evolution of the subhalo masses in Figure^ 
we can identify the infall times of the subhalos into the Milky Way 
halo. In CDM, the subhalos that match the Fornax kinematics and 
luminosity have infall redshifts of z = 1.25 and 2 = 1.37, cor¬ 
responding to lookback times of ~ 9 Gyr. For the 2.3 keV WDM 
model the infall redshift is 2 = 1.02, corresponding to a lookback 
time of ~ 8 Gyr. For DDM model the infall redshift is at 2 = 0.28, 
corresponding to a lookback time of ~ 3 Gyr. So the DDM host 
subhalo is predicted to fall into the Milky Way more recently than 
the CDM and WDM candidates. 

It is informative to compare the lookback times that we de¬ 
duce to_greyiot^estimates of the infall times of Milky Way satel¬ 
lites. [Rocha^T^lloT^) estimate infall times using subhalo Galac- 
tocentric positions and orbital motions from Via Lactea II simula¬ 
tion. Using these criteria Fornax is most likely to have fallen into 
the Milky Way halo ~ 5 — 9 Gyr ago. So even though our crite¬ 
ria for identifying Fornax subhalo cand i dates in CDM simulations 
are different from those in iRocha et alj ( 1201 2h . the infall times are 
in good agreement. We note that the Aquarius simulations that we 
utilize and our stricter matching criteria using luminosity and stel¬ 
lar kinematic data allows for a lar ger range of merger histories than 
considered in iRocha et alj ( l2012h . who utilize a single Milky Way 
halo. 


4.3.3 Star-Formation Quenching Mechanism 

The quenching of star formation in Milky Way dSphs may be re¬ 
lated to their infall times. For example ram-pressure stripping can 


remove the cold gas at the center of the satellites as a result of 
the high-speed inter action with the hot gas halo of the host halo 
dGunn & Gotlj|l97^ . Ram pressure stripping has been invoked as 
the quenching mechanism for Milky Way and M31 dwarf satellites 
with Mi, 10 ^ Mp) with an ex tremely short quenching timescale 
~ 2 Gyr dFlllingham et aljEoiSh . 

As shown in Figure[8l Fornax is observed to have an enhanced 
star formation activity ~ 3 -4 Gyr ago and quenching ~ 2 Gyr ago 
dColeman & de Joir^boOSl) . If we match these timescales with the 
infall times that we determined above, the enhanced star forma¬ 
tion event is most consistent with the infall time of our DDM sub¬ 
halo candidate. The g rey shaded area sho ws the possible Fornax 
Infall time range from iRocha et alj d2012h . This covers a large in¬ 
fall time range including those predicted by the CDM and WDM 
simulations. However, it is interesting to note that other recent For¬ 
nax st ar formation history studies either indicate similar peak at ~ 
4 Gyr dde Boer et al.l201% o r draw different conc lusions with a star 
formation peak at ~ 8 Gyr ddel Pino et al.ll20i^ . Deep photomet¬ 
ric studies of Fornax will help resolve this discrepancy and provide 
more insight on the connection of the Fornax star formation history 
with its infall time. 

To explain the Fornax star formation history in the context of 
the earlier subhalo infall times predicted in our CDM and WDM 
models, we can consider a scenario in which the Fornax star for¬ 
mation was quenched when it merged into Milky Way halo ~ 9 
Gyr ago. A close encounter with the Milky Way at the last peri- 
galactic passage would then trigger star formation. This mecha- 
n ism has been invok ed to explain the star formation history of Leo 
I dMateo et alj2Q08h . which had a burst in its star formation ~ 3 — 4 
Gyr ago. Indeed proper motion and line-of-sight velocity measure¬ 
ments indic ate that Leo I is l ikely on a fairly eccentric, nearly un¬ 
bound orbit dSohn et alj2013h . In contrast with Leo I, Fornax likel y 
has a much less eccentric orbit dLux et al.l2010l : IPiatek et al.l2007h . 
so a much closer encounter with the Milky Way may be more dif¬ 
ficult to Invoke in the case of Fornax. More precise measurement 
of orbital velocity and detailed orbit reconstruction will shed more 
light on this issue. 
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5 CONCLUSIONS AND DISCUSSION 

We have fit the kinematics and photometry of the Fornax dSph to 
subhalo potentials predicted by numerical simulations of different 
dark matter models. CDM subhalo candidates are typically in the 
mass range of 10* — 10® Mq. WDM and DDM models predict on 
average larger mass subhalo hosts of Fornax. 

The diversities of subhalo candidate properties in different 
dark matter models and the different time evolution of dark matter 
free-streaming effects lead to diversities in subhalo formation his¬ 
tory, and also likely their stellar formation history. With this in mind 
we utilize simulation merger trees to investigate the formation his¬ 
tory of the subhalos that match the Fornax kinematics. We imple¬ 
ment simple models to match the Fornax luminosity to dark matter 
halo mass in order to understand possible star formation histories 
of these Fornax subhalo candidates. Under the assumption that the 
mass of the most massive progenitor is correlated with the current 
stellar content in each galaxy, we derive current stellar masses for 
each subhalo. 

For CDM subhalos, our best candidates require the ratio be¬ 
tween the peak mass and current mass to be ^ 10. This suggests 
that these systems experience significant tidal stripping after they 
fall into Galactic halo and are subject to substantial dark matter 
mass loss. It also implies that their infall times are ^ 9 Gyr ago, so 
that they fall well inside the galactic halo where the tidal stripping is 
ef ficient. Interestingly, similar findings also have been pointed out 
bv ICooper et al.l j2010l) using Aquarius simulations with a “particle 
tagging” technique to simulate Galactic stellar halos. Using Fornax 
and Carina as two examples, these authors find that they can only 
match their observed dwarf galaxy surface brightness and veloc¬ 
ity dispersion profiles simultaneously by choosing model satellites 
that have suffered substantial tidal stripping. 

Our WDM subhalo candidates show only mild deviation from 
their CDM counterparts, in that they have similar peak mass to cur¬ 
rent mass ratio and infall time. A caveat to this analysis is the real 
stellar content might be different, and full hydrodynamic simula¬ 
tions are needed to calibrate galaxy formation process in WDM. 
For the DDM candidates, using similar assumptions our subhalo 
candidates have relatively small peak mass/current mass ratio (~ 
3) and they are more massive at z = 0 than the CDM candi¬ 
dates. Therefore the DDM candidate only experiences mild mass 
loss since merging into the Galactic halo approximately 3 Gyr ago. 

We note that our results are not yet able predict the existence 
of visible stellar tidal tails. Since galaxies are embedded deep in 
the center of the halo potential, most of the tidal stripping effects 
are expected to be on the DM component. Determining whether 
or not sta rs are stripped out re quires N-body simulations with star 
particles. ICooner et alj ( l2010h have shown that 20% of the stellar 
mass and 2% of DM mass of thei r Fornax candidate re mains at 
present time. Also recent work from iBattaglia et aTI (1201 Sll have ex¬ 
amined the tidal effects on Fornax along its possible orbits. Current 
and forthcoming wide field deep imaging surveys may be sensi¬ 
tive to faint ste llar tidal tails and can provide an observational test 
of our models. [Bate et aP ( 1201 Sll use data from VLT Survey Tele¬ 
scope (VST) ATLAS Survey to study a region of 25 square degrees 
centered on Fornax. They have excluded a shell structure outside 
of Fornax’s tidal rad ius that was reported from a previous study 
dColeman et al.ll200^ . In the near future the Dark Energy Survey 
(DES) will cover a wide area around Fornax with much improved 
imaging depth, and possibly provide a test of the tidal stripping 
hypothesis. 

We find that most of the subhalo candidates are subject to the 



Lookback ti(me [Gyr] 

Figu re 8. Qualitative Fornax s tar-formation histories (black solid lines, 
from| Coleman & de Jond d2008h l compared to estimated infall times from 
iRocha etakfllOl^ using Galactocentric positions and orbital motions (grey 
filled area) and our results (vertical dashed lines). The two navy dashed lines 
indicate the infall time for the two CDM candidates from the Aquarius sim¬ 
ulations that match both the Fornax luminosity and stellar kinematics. The 
aqua line is for the WDM 2.3 keV candidate, and the orange line is for 
the DDM can didate. Notice the D DM candidate infall time is significantly 
lower than the lRocha et alj l2012l) prediction and our CDM and WDM can¬ 
didates, and also it happens at the time when the Fornax star formation rate 
starts to drop after a burst 3-4 Gyr ago. All others predict infall times before 
the enhanced star formation event. 

cosmic reionization UV-background so that their star formation is 
highly suppressed. In our sample more than half of the CDM candi¬ 
dates are likely not visible at all due to reionization. For the WDM 
model the reionization history is more uncertain because high- 
redshift structure formation is significantly delayed. However, for 
DDM the early structure formation follows CDM until late times 
when the decay process becomes significant, so a similar reioniza¬ 
tion history can be applied. Most of our DDM candidates are not 
affected by the UV-background since they are already massive at 
z=6. This is largely because for the same Fornax kinematic data, 
DDM candidates tend to be more massive because they have less 
concentrated DM profiles than their CDM counterparts. 

The difference in Fornax infall times in the different models 
may imply a different star formation quenching mechanism if en¬ 
vironmental effects play an important role in dSph formation pro¬ 
cesses. From the observed Fornax star formation history, there is 
an enhanced s tar formation event at 3-4 Gyr ago and it is quenched 
~ 2 Gyr ago dColeman & de Jon3l2008l) . For those subhalo candi¬ 
dates that have infall times well before the star formation peak, star 
formation may be triggered by a close passage to the Milky Way. 
For those that have infall times after the star formation peak, the in¬ 
fall into Galactic halo may cause rapid gas loss due to ram-pressure 
stripping. More precise predictions can be made by detailed Fornax 
orbital motion reconstruction from its proper motion and line-of- 
sight velocity measurement, which may be possible with the GAIA 
satellite. Larger sets of simulations that sample more possible halo 
formation histories in different dark matter scenarios will also help 
to confirm our results. 
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APPENDIX A: RESOLUTION TESTS AND 
CORRECTIONS 

As pointed out in iBovlan-Kolchin et al.l ( 1201 2h where they used 
Aquarius simulations to compare dynamic properties of simulated 
galactic subhalos with observed Milky Way satellites, the effect of 
force softening will introduce errors to the subhalo potentials in the 
central region. The reason is that the force softening scales of our 
simulations range from e = 65.8 - 72.0 pc. The effect of force soft¬ 
ening will reduce the density on scales of ~ 3 e. However we need 
to work with a cumulative mass profile that will carry the effect to 
larger scale. In order to estimate and further correct for this effect, 
here we utilize the Aquarius A level 1 simulation, which has the 
same halo realization as the Aquarius A level 2 simulation but with 
much smaller force softening length (e = 20 pc) and particle mass 
(mp = 1.712 xlO® Af©), to test different correction methods to 
reconstruct the mass distribution at subhalo center. 

The correction method adopted here is similar to the one in 


iBovlan-Kolchin et al.l ( l2012h . We have tried out a few different fit¬ 
ting ranges and correction ranges to find the one that gives the best 
results. The procedure is as follows. We fit the subhalo density pro¬ 
files with an Einasto profile, which has the following form: 


p{r) = p -2 exp 



a 



(Al) 


The fitting range is [rinner ,rupper], where rupper is the smaller 
of 3 kpc or 1.5 Tmax and we try a few different values for Tinner- 
We then replace the mass within Tinner with the accumulated mass 
generated from the Einasto profile fit: 


A,«,)^ 4rrr-,p-, exp( 


3 Ina + 2 — ln8 \ 

-1-h 


3 2( 
a a\T-2/ 


(A2) 


and use the mass profile from simulations outside of the Tinner- 

In Figure IaTI we compare the Fornax data fitting results us¬ 
ing different corrections with ditferent Tinner on Aquarius level 
2 subhalos to their Aquarius level 1 counterparts. Here we show 
three different Tinner choices: no correction (Tinner = oo, red his- 
to grams in FigurelAlll. rinner = 300 pc (close to choice of 291 pc 
in Bovlan-Kolchin et alj(l2012h . blue histograms in Figure lATt . and 
200 pc (green histograms in Figure lATt . We take the ratio of the fit¬ 
ting parameter value from Aquarius level 2 subhalos with their level 
1 counterparts. So the ideal case will be they lie very close to 1, in¬ 
dicating that the corrections faithfully recover the mass distribution 
in higher resolution simulation runs. 

We can see that some of the parameters are not affected much 
by the correction. Examples are c and ro, since they are related to 
stellar distribution in the outer region (for example, ro ~ 1 kpc) 
and are not sensitive to the correction in the inner region (^ 300 
pc). The of ths photometry data fit also is not affected signif¬ 
icantly, since most of the photometry data points are distributed 
outside ^ 300 pc. The worst of all parameters is a, and the amount 
of deviation seems to be large no matter which correction we ap¬ 
ply. The reason is that the data do not provide good constraints on 
the parameter a. The uncertainty of the MCMC result for a is usu¬ 
ally very large and highly correlated with other parameters such as 
the parameter that determines the normalization of the photometry 
data. So the wide distribution of the parameter a for different cor¬ 
rections is likely not due to the way the dark matter potential is cor¬ 
rected, but due to the large uncertainties from the poor constraints 
and parameter correlations. From the overall behavior of these three 
different correction methods, we can see that the green histograms, 
for which the Einasto profile correction is applied within 200 pc, 
show the best agreement with Aquarius level 1 counterparts. 

In Figure lAZl we show the fractional difference in circular ve¬ 
locity between Aquarius level 1 subhalos and their level 2 counter¬ 
parts after corrections at two radii: r=200 and 300 pc. At 300 pc 
the difference is within 20% and at 200 pc the difference is within 
30% for all correction methods. However, before applying any cor¬ 
rection, the difference is negative (red histograms) and for the cor¬ 
rection applied within 300 pc (blue histograms) the difference is 
biased toward positive values. The green histogram, which is for 
the correction within 200 pc, is peaked at zero and with difference 
within 10-15%. This indicates that before any corrections the cen¬ 
tral dark matter potentials (Vcir = \/GM{r <)/t) are reduced 
due to softening scale effects for Aquarius level 2 simulations, and 
the suppression is up to abo ut 20%. However, if w e apply the cor¬ 
rection method suggested in IBovlan-Kolchin et alj (120121). we can 
’’over-correct” the potentials up to 20-30% at 200-300 pc. There¬ 
fore our tests show that a mild correction within 200 pc with the 
fitted Einasto profile mass will generate the best results for remov- 
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Figure Al. The distribution of the fitting results for Aquaiius A-2 simulations with different softening scale correction methods : no comection (red histogram), 
blue histogram with rinner= 300 pc, and green histogram rinner='^00 pc. Here Tinner is within what radius a Einasto profile fit is used for accumulated 
mass profile. The x-axes show the ratio of fits to their Aquarius level 1 counterparts. 
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Figure A2. The distribution of circular velocity (Vcir) fractional difference 
for different Aquarius A resolution simulations and different softening scale 
correction methods. The x-axis shows the fractional difference of the Veir- 
The color scheme is the same as Figure I aTI 


ing the softening scale effects. We thus adopt the correction all the 
simulations in this study. 


APPENDIX B: EFFECTS OF WMAP 1 & WMAP 7 
COSMOLOGICAL PARAMETERS 

Here we test the effects of differe nt cosmological parameter s 
on our study. Previous stu dies (e.g. IZentner & Bullocfl (l2003h : 
[Polisenskv & Ricottil ( |2014|) ) have shown that different cosmolog¬ 
ical parameter values, especially as and ris, have non-negligible 
effect to the substructure density profiles. The simulation sets used 
in our study have adopted different cosmological parameter val¬ 
ues. We compare the Fornax Jeans equation fitting results from the 
Aquarius A Level-2 simulation that is re-simulated with a WMAP 
7 cosmology (Aq-A2-w7) with those from the original Aquarius A 
Level-2 simulation that is simulated with a WMAP 7 cosmology 
(Aq-A2). One of the most significant differences of their choice of 
parameters lies in the difference in cts (cts = 0.9 for Aq-A2 while 
(Ts = 0.81 for Aq-A2-w7). In the following we will compare the 


fitting results and discuss the impacts due to the differences in cos¬ 
mological parameters. 

In Figure ISTI we show the present-time subhalo properties of 
the Fornax candidates and their circular velocity curves. Those are 
candidates that fit Fornax kinematics but not necessarily fit the For¬ 
nax luminosity. In both simulations we find a comparable number 
of good-fits. Here we only show those we can find counterparts in 
both simulations. Their maximum circular velocities (Vmax) and 
subhalo mass (Msub) span very similar range. However, it is shown 
in the right panel in Figure I bT] that those subhalos from Aq-A2 are 
slightly more co ncentrated at the center than t hose from Aq-A2-w7. 
As discussed in IPolisenskv & Ricottil ( l2014h . halos in simulations 
with higher ag form earlier and thus are more concentrated. How¬ 
ever, we can see from the distribution of candidate subhalo proper¬ 
ties that the changes due to the differences in WMAP 1 and WMAP 
7 cosmology are much smaller than the difference due to different 
dark matter properties. When we compare Figure lSTI with Figure|4] 
we can see that WDM and DDM simulations generate subhalos that 
exhibit very different properties that occupy a completely different 
subhalo property region. It is thus clear that high-quality Fornax 
data have provided good constraints on subhalo present-time prop¬ 
erties. 

We also investigate the differences in halo formation history. 
In Figure IB^ we check the distribution of Mpeak and mass at z=6 
of Fornax candidates from these two simulations. We find that the 
difference is mild and does not affect our conclusion. For exam¬ 
ple, for both simulations more than half of the candidates may have 
suppressed star-formation due to the UV-background from reion¬ 
ization. Also for both simulations only one candidate matches the 
Fornax luminosity while others are much dimmer. Thus we con¬ 
clude that the difference in WMAP 1 and WMAP 7 is mild and it 
has limited effects on both present-time properties and halo forma¬ 
tion history. The differences due to dark matter properties, which is 
the focus of our study, exhibit much stronger effects. 
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Figure Bl. Subhalo properties at z=0 : maximum circular velocity (Vmax, left panel), and subhalo mass {Mgub^ middle panel), and circular velocity profiles 
(right panel) for the Fornax good-fits from AqA WMAP 1 (red dash-dotted lines) & AqA WMAP 7 (blue solid lines) simulations. Here we only show those 
we can find counterparts in both simulations. 




Figure B2. The Mpeak and mass at z=6 {Mz=6) distribution of the Fornax 
good-fits from the AqA WMAP 1 (blue solid histogram) and AqA WMAP 
7 (red dash-dotted histogram) simulations. The shaded area in the left panel 
indicates the progenitor mass range with baryon content that is subject to the 
UV background suppression. The yellow shaded area in the right panel indi¬ 
cates the region where the luminosity of these objects are co nsistent with the 
Fornax luminosity using abundance matching methods from lBehroozi et ^ 
<2013l) . 














































































